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SUMMARY 

Classical undamped acoustic-wave theory was used to determine analytical re- 
lations among sinusoidal perturbations of pressure and flow at the ends of a 
hydraulic- transmission line having a closed-end branch of arbitrary length a - 
tached at an arbitrary point. Experimental data were obtained for the equila er- 
ai case (a branch half as long as the main line connected to the main line at the 
midpoint) at mean flow speeds of 5 to 10 feet per second. Measured pressure- 
perturbation ratios agreed closely with analytical predictions. At frequencies 
for which the branch length was an odd multiple of 1/4 wavelength, waves m the 
main line were almost completely reflected at the junction point. 


INTRODUCTION 


Acoustic pressure and flow perturbations in long hydraulic lines are often 
of importance in the operation of various fluid systems. Among these are hydrau- 
lic control systems (ref. l) and rocket propellant systems (ref. 2). Because of 
the need to predict acoustic effects in fluid systems for rocket and space appli- 
cations, the NASA Lewis Research Center has undertaken a study of acoustic dis- 
turbances in flowing-liquid systems. An investigation of acoustic waves in a 
long straight line rigidly supported at one end is reported in reference 3. e 
effect of longitudinal pipe motion produced by unbalanced pressure forces on an 
elbow is treated in reference 4. 


This report presents the results of a study concerning the effects of a 
closed branch on acoustic disturbances in a line. A general analysis based on 
classical acoustic-wave theory is presented. This analysis shows the effect of 
a closed branch of arbitrary length attached to the main line at an. arbitrary 
point when the line is terminated in an arbitrary impedance. ^Experimental re- 
sults are given for a particular configuration in which a 34-loot closed branci 
was attached at the midpoint of a line G8 feet long that terminated in a resis- 
tive impedance. The main line and branch were of stainless-steel tube 1.00 in i 
in outside diameter with a 0. 065-inch wall. The experiment covered a range of 
disturbance frequencies from 0.5 to 02.0 cycles per second, mean flow speed in 
the main line from 5 to 10 feet per second, and mean line pressure iron 50 to 
225 pounds per square inch gage. The average amplitude of the sinusoidal pertur- 
bation was about 5 pounds per square inch for pressure and 0.13 foot per second 


for fluid velocity. The flow of the fluid (JP-4 fuel) was modulated by means of 
a hydraulic servovalve at the upstream end and was restricted by an orifice at 
the downstream end. 


EXPERIMEIMTAL SYSTEM 
Apparatus 

Flow. system . - The essential parts of the open-loop pumped- return flow sys- 
tem used in the experiment are shown in figure 1. The fluid (JP-4 fuel) was 
forced through the test line by a gear pump, and the mean flow rate was measured 
by the use of a rotameter. Hydraulic accumulators were placed between the pump 
and the test line to provide steady supply pressure. The discharge from the test 
line was submerged in a constant-height vented tank. Fluid was returned to the 
supply tank by intermittent operation of the return pump. 

Flow disturbance generator . - Sinusoidal perturbations of flow and pressure 
were induced in the system by means of an electrohydraulic throttle servovalve 
located just upstream of the test line. The throttle was oscillated sinusoidally 
about a partly open mean position in response to an alternating voltage. A 
porous -metal filter was placed between the throttle and the test line to reduce 
turbulence. 

Test line . - The test line, made of stainless- steel tubing having a 
1.000-inch outside diameter and a 0.065-inch wall thickness, was fitted with a 
right-angle tee at the midpoint. Attached to the tee was a 34-foot branch line 
of the same tubing that was closed at the end by a rigid cap. Provision was made 
for bleeding the line to get rid of air. The closed end of the branch and both 
ends of the main line were clamped to rigid supports in order to avoid pipe- 
motion effects of the type described in reference 4. The downstream end was ter- 
minated in an orifice plate containing 34 holes 0.040 inch in diameter. 

Instrumentation . - Instrument sensing elements to measure pressure perturba- 
tions were located at stations A, B, C, D, E, and F, as shown in figure 1. The 
pressure sensors were commercial flush- diaphram units. Pressure perturbations in 
pounds per square foot were obtained by means of static calibration of each unit 
with its associated amplifier. The in-phase and quadrature components of the 
pressure perturbations were read from a commercial transfer- function analyzer, as 
in the investigations previously reported (refs. 3 and 4). 


Procedure 


The flow system was operated at constant mean flow rate and pressure until 
conditions had stabilized. The branch line was bled at the closed end to ensure 
that no air was trapped in the system. Mean flow rate and pressures at stations 
A and E were recorded. The throttle was operated at frequencies of 0.5 to 82.0 
cycles per second. The amplitude of the throttle area variation was maintained 
constant over the entire range of frequencies for each run and was kept small 
relative to the mean open area to avoid nonlinear effects. At each frequency, 
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the values of the in-phase and quadrature components of the several pressure 
perturbations were read on the transfer-function analyzer. The mean pressure 
drop across the downstream orifice was constant during each run, but was changed 
from run to run to vary the orifice impedance. The exit orifice impedance was 
determined for each run from the slope of the steady-state pressure-flow curve. 

The following dimensionless pressure-perturbation ratios were computed from 
the data for each frequency: 

(!) P E /P A 

(2) Pp/P A 

(3) Pb/P a 

(Symbols are defined in appendix A. ) 


ANALYSIS 

The system to be analyzed consists of a straight main pipe of length Z and 
a branch of the same diameter and of length kZ, which is attached to the main 
line at a distance xZ from the downstream end, as shown in figure 2. The down- 
stream end of the main line E (fig. 2) is terminated in an arbitrary load im- 
pedance Z E . The end of the branch F is closed by a rigid cap. The positive 
direction of flow is chosen from A to E in the main line and from D to F 
in the branch. With this convention for the direction of the exit flow perturba- 
tion Qg, the exit load impedance is defined by 



This equation represents the impedance looking into the terminating structure 
from the end of the pipe. 

In the following analysis, dimensionless complex numbers are used to repre- 
sent pressure and flow perturbations and their ratios. The amplitude and phase 
of the various sinusoidal perturbations and their ratios are related to the cor- 
responding complex numbers according to the usual conventions. Dimensionless 
pressure and flow perturbations and impedance are defined as follows: 


p 

p - — 2 

(2) 



q-4 

(3) 

H cS 


_ P _ zs 

(4) 

q ~ pc 
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A dimensionless frequency parameter defined by 


( 5 ) 

represents the phase shift in the distance l for a single sinusoidal. wave train 
of frequency f and propagation speed c. 


It is assumed that undamped sinusoidal acoustic waves exist in the system, 
that the pipe does not move longitudinally in response to the pressure perturba- 
and that the mean fluid flow speed in the main line is negligible compared 
to the sonic speed. The acoustic line equations (derived in ref. 4) relating the 
pressure and flow perturbations at the two ends of a straight section of uniform 


pipe are, for section CE (fig. 2), 

P C = cosh ix P + d E sinh ix P (6a) 

dp = d E cosh ixp + p E sinh ixf3 (6b) 

PE = Pc cosk i x P - Op sinh ixp (6c) 

d£ = dp cosk ixP - Pq sinh ixp (Sd) 

The following similar equations apply for sections AB and DP, the phase 
shifts being (l - x)p and kp, respectively: 

PA = PB cosh i(l " x) p + qg sinh i(l - x)p (7a) 

q A = q B cosh i(l - x)p + p B sinh i(l - x)p (7b) 

p B = P A cosh i(l - x)p - q A sinh i(l - x)p (7c) 

d B = d A cosl1 i(l - x)p - p A sinh i(l - x)p (7d.) 

Pp = Pp cosh ik P + d F sinh ik P (8a) 

dp dp cosh ikp + Pp sinh iicp (Bb) 

P F = p D cosh ikp - q D sinh ikp (e c ) 

dp = dp cosh ikp - Pp sinh ikp (0<j) 

The impedance relation at the downstream end (station E) is 


p E 

^ “ ZE (9) 
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The impedance at the closed end of the branch (station F) is infinite hence , 




( 10 ) 


Continuity of flow and pressure at the junction point is expressed by 

= ^ + q D 


and 


PB = Pc " Pd 

The corresponding impedance relation is 


1_ 

Z B 


Z C Z D 


( 11 ) 


( 12 ) 


(13) 


IS 


where z-p is the impedance looking into the branch from the junction, z c 
looking downstream from the junction in the main line, and z B is looking into 
the junction from the upstream side. 

Transfer functions relating the various perturbations were obtained from 
equations (6) to (13) (mathematical details in appendix B). The relations among 
the perturbations at the two ends of the main line are conveniently expressed in 


LS 

of p, 

, k. 

and 

X 

defined as follows 

T 

= cos 

P - 

tan 

kp 

cos 

xp 

sin(l - 

x)0 

U 

- sin 

P - 

tan 

kp 

sin 

xp 

sin(l - 

X)0 

V 

~ cos 

P - 

tan 

kp 

sin 

xp 

cos(l - 

x)p 

w 

ss sin 

0 + 

tan 

kp 

cos 

xp 

cos(l - 

x)p . 


(14) 


In terms of these functions, the entrance impedance (looking downstream at 
point A) is 


Tz-g + iU 
z -^ _ V + iz-gW 


The interterminal pressure ratio is 




•E 




(15) 


( 16 ) 


The interterminal flow ratio is 
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Other ratios of possible interest are 


and 


pa 


Tz-g + iU 


P E _ Z E 
<1^ V + iz^W 


(18) 


(19) 


Equations (14=) to (19) define the dynamic behavior of the system in terms of the 
pressure and flow perturbations at the terminals of the main line. In the ex- 
perimental investigation, pressure perturbations were also measured at the junc- 
tion point and at the closed end of the branch. Relations (derived in appen- 
dix B) involving Pg and Pp required for comparison with experimental results 
are as follows : 


Pp 

z-g cog x(3 -f I sin x(B 

(20) 

Pa = 

cos kp(TZg + iU) 

P B 

z E cos xJ3 + i sin x|3 

(21) 

pa " 

Tzg + iU 


RESULTS AND DISCUSSION 

Pressure perturbations measured at stations B, C, and D on the three sides 
of the tee agreed within the limits of experimental error and, thus, verified the 
assumption of identical pressure perturbations in all branches at the junction 
point. Calculated and experimental pressure perturbation ratios are presented in 
figures 3 to 5. The curves were calculated from equations (IS), (20), and (21); 
data points were determined from experimental measurements as described in the 
section Procedure. The ratio of downstream to upstream pressure is shown in fig- 
ure 3. Figure 4 gives the ratio of pressure at the closed end to upstream pres- 
sure, and figure 5 gives the ratio of pressure at the junction point to upstream 
pressure. The experimental and analytical results are in agreement, which shows 
that classical undamped acoustic-wave theory provides a good description of the 
actual fluid wave motion for the conditions of the experiment (low- viscosity 
liquid and mean flow speed small compared to sonic speed). 

The expressions for impedance and flow- perturbation ratios given in equa- 
tions (15), (17), (18), and (19) could not be checked experimentally because flow 
perturbations were not measured, but the successful verification of equa- 
tions (16), (20), and (21) is reason for confidence in the applicability of the 
others. Values of the upstream impedance calculated from equation (15) for the 


experimental configuration (k = x — 0.5) are shown in figure 6. 

Because of the complete reflection of waves at the closed end of the branch 
line, the system exhibits strong resonances at frequencies making the branch 
length an odd multiple of l/4 wavelength (k0 = */2, 3rt/2, etc.). At these fre- 
quencies the impedance looking into the branch is practically zero, and the pres- 
sure perturbation at the junction is, therefore, infinitesimal. Under these con- 
ditions, the branch accepts practically the entire flow perturbation entering the 
tee from the upstream line, and virtually no wave motion exists in the main line 
downstream of the junction. 


CONCLUDING REMARKS 

The results of this experiment and the two preceding investigations at this 
research center (refs. 3 and 4) indicate that undamped acoustic-wave theory de- 
rived for stationary fluids satisfactorily describes disturbance propagation in 
liquids flowing within long propellant lines, provided that any motion of the 
line itself is taken into account. Three different configurations, a straight 
line, a 90° L, and an equilateral T, were tested. 

Similar analytical methods based on undamped wave theory would be expected 
to apply to small perturbations in other propellant-line configurations under 
the following conditions : 

(1) Mean flow speeds up to about 20 feet per second 

(2) Length-to-diameter ratios to about 1000 

(3) Frequencies to about 100 cycles per second 

(4) Fluids of low kinematic viscosity and high sonic speed such as water, 

kerosene, alcohol, liquid metals, and subcooled liquefied gases 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, May 10, 1963 
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APPENDIX A 


SYMBOLS 


c 

f 

k 

l 

P 

P 

Q 

q 

s 

T,U,V,W 

x 

Z 

z 

0 

P 

Subscripts : 


sonic speed, ft/sec 

frequency of sinusoidal perturbation, cps 

ratio of branch to main-line lengthy dimensionless 

length of main line, ft 

perturbation pressure (complex number) , lb/sq ft 

dimensionless perturbation pressure (complex number) , P/pc^ 

perturbation flow (complex number), cu ft/sec 

dimensionless perturbation flow (complex number) , Q/cS 

inside area of cross section of pipe, sq ft 

function of p, x, and k defined by eq. (14), dimensionless 

position coordinate of junction (ratio of distance from exit to 
junction to main-line length), dimensionless 

acoustic impedance, (lb) (sec) /ft 5 

dimensionless acoustic impedance (complex number), ZS/pc 
frequency parameter, 2*fZ/c, dimensionless 
mean fluid density, slugs/cu ft 


A,B,C,D,EjF reference stations located as shown in fig. 2 
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APPENDIX B 


DERIVATIONS 

Transfer functions relating the various perturbations of interest may be 
tained conveniently by determining the ratio of each to the exit pressure per 
tubation pg. 

From equations ( 6 ) , ( 9 ) , and ( 12 ) 

2= = 5 = 2 = cosh ixp + slnb ±x & 

P E % P E Z E 


and 


_C _ cosh ixfi + sinh ixp 


■E 


J E 


From equations ( 8 ) and (10) 


Pt) 

— = cosh lkp 


and 


*0 


sinh ikp 


The ratios Pp/pg and dp/pp are found as follows: 

, . n . sinh ixp 

cosh xxp + 

Pp Pf Pd Z E 


P E P D P E 


2 = -X - tanh ikp (cosh ixp + 


cosh ikp 


-E 


Pg P E 




\ 


sinh ixp 
Z E / 


From, equation (ll) 


+ ^D 

P E P E P E 


Substituting from equations (B2) and (B6) gives 


ob- 

(Bl) 

(B2) 

(B3) 

(B4) 

(B5) 

(B6) 

(B7) 
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qg i 

— = sinh ixp + tanh ikp cosh ixp + — (cosh ixp + tanh ikp sinh ixp) 


When equation (7a) is divided, by j>g, 

P a P-n q^n 

— = — cosh i(l - x)p H sinh i(l - x)p 

P E P E P E 

Substituting from equations (Bl) and (B8) and simplifying result in 


— = cosh ip + tanh ikp cosh ixp sinh i(l - x)p 


+ — ["sinh ip + tanh ikp sinh ixp sinh i(l - x)pl (BIO) 

Z E L J 

Replacing hyperbolic with trigonometric functions and using the notation of equa- 
tions (14) give 


IA = T + iU 


When equation (7b) is divided by p^. 


q a qg p-n 

— = — cosh i(l - x)p H sinh i(l - x)p 

P E P E PE 


(Bll) 


Substituting from equations (Bl) and (B8), simplifying, and using the func- 
tions defined by equations (14) give 


jA _ _V_ 
P E Z E 


Equation (19) is the reciprocal of (B12). 


(B12) 


q^ v + i zgW 


le product of equations (16) and (19) is 


Tzj. + iU 

s a = vTTZw 

E 
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Multiplying equation (16) by equation (9) yields 


Eh = Tzxp + iU ( 18 ) 

*E E 

The product of equations (9) and (B12) is 

= V + iZpW (17) 

Dividing equation (B5) by equation (16) and replacing hyperbolic by trigonometric 
functions give 

Pp z g cos xp + i sin xp (20) 

~ cos kp(Tz E + iU) 


The ratio p^/p^, obtained in the same way from equations (Bl) and (16)., is given 
by 


p^ z_ cos xp + i sin xp 
Jd iii 

Pa = Tz e + 1U 


( 21 ) 
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